Magnetic properties of single atoms of Fe and Co on Ir(lll) and Pt(lll) 
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In using the fully relativistic versions of the Embedded Cluster and Screened Korringa-Kohn- 
Rostoker methods for semi-infinite systems the magnetic properties of single adatoms of Fe and Co 
on Ir(lll) and Pt(lll) are studied. It is found that for Pt(lll) Fe and Co adatoms are strongly 
perpendicularly oriented, while on Ir(lll) the orientation of the magnetization is only out-of-plane 
for a Co adatom, for an Fe adatom it is in-plane. For comparison also the so-called band energy 
parts of the anisotropy energy of a single layer of Fe and Co on these two substrates are shown. The 
obtained results are compared to recent experimental studies using e.g. the spin-polarized STM 
technique. 

PACS numbers: 73.20. At, 72.10.Fk, 73.22.-f, 75.30.Hx, 73.20.Hb 



The potential application in non-volatile data storage 
devices is one of the driving forces behind research into 
magnetic nanostructures. In state-of-the-art hard disk 
drives a collection of a few hundred of single-domain par- 
ticles (grains) are used to hold one bit of information. If 
materials can be manufactured which exhibit sufficiently 
large anisotropics and thermal stabilities it may become 
possible to store one bit in a single grain [1 . Such storage 
devices will require magnetic structures of precise atomic 
arrangement, as - if in addition the lateral dimensions of 
grains are further reduced - the influence of the perime- 
ter atoms becomes increasingly important [2j |3] and, as 
is known, from previous studies the magnetic properties 
of each atom in a nanostructure are highly influenced by 
its local environment [TJ |2j |3j [4] . 

Using Scanning Tunneling Microscopy structures can 
be precisely tailored and their magnetic properties de- 
termined. In recent Scanning Tunneling Spectroscopy 
experiments [31 |5] it has become possible to measure 
not only the Lande g-factor of individual atoms but also 
their magnetic anisotropy. The findings suggest that the 
anisotropy energy of a single atom may eventually be 
large enough to use the magnetic state of an atom as 
a storage unit, pushing the ultimate limit for data stor- 
age density even further. Since the magnetic properties 
of small clusters and single adatoms differ strongly from 
those of bulk systems and even monolayers - e.g. show- 
ing a much enhanced magnetic anisotropy energy - they 
do not only generate interest for their technological rele- 
vance but also from a fundamental point of view. 

In this paper we present a study of the magnetic mo- 
ments and the angular dependent band energy part of the 
magnetic anisotropy energy of single atoms of Fe and Co, 
which - in order to investigate the influence of different 
substrates - have been deposited on Pt(lll) and Ir(lll). 
The calculations have been performed by means of the 
Embedded Cluster Method (ECM), a scheme based on 
the fully relativistic Screened Korringa-Kohn-Rostoker 



(SKKR) method, in which we can treat impurities em- 
bedded into a two-dimensional translationally invariant 
semi-infinite host. This approach makes use of multiple 
scattering theory in which the electronic structure of a 
cluster of embedded atoms is described by the so-called 
scattering path operator given by the following Dyson 
equation [6j[7]: 

r c (e) = r h (e)[l - (t^(e) - (e))^)]" 1 (1) 

where r c (e) and r^(e) account for all the scattering events 
within the embedded cluster and the host, t~ l and t^ 1 
denote the single-site scattering matrices for the "impu- 
rity" and for the host atoms, respectively. Once r c (e) 
is known all corresponding local quantities, i.e., charge 
and magnetization densities, spin and orbital moments, 
as well as the total energy can be calculated. 

In order to perform self-consistent calculations within 
local density functional theory [8], for the calculation of 
the ^-matrices and for the multipole expansion of the 
charge densities (needed to evaluate the Madelung po- 
tentials), a cutoff for the angular momentum expansion 
of /max = 2 was used. The potentials were treated within 
the atomic sphere approximation (ASA). In all self con- 
sistent calculations the orientation of the magnetization 
was chosen to point uniformly along the surface normal (z 
axis) . Structural relaxations of the cluster-substrate dis- 
tance, which may in principle affect the magnetic prop- 
erties [9j [TUJ [TTJ H2] , have been neglected. The host 
and the impurity sites refer to the positions of an ideal 
fee lattice with the experimental lattice constants of Pt 
(a = 3.92 A) and Ir (a = 3.84 A). The self-consistent 
calculations were performed using 102 fen -points in the 
irreducible Surface Brillouin zone integrations and 16 en- 
ergy points for the energy integrations along a semicir- 
cular contour in the complex energy plane by means of 
a Gaussian quadrature. To guarantee that all the per- 
turbed host atoms with non-negligible influence on the 
calculated properties are taken properly into account, we 



2 




I i I i I i L 

12 3 
shells 



FIG. 1: Calculated spin (top) and orbital (center) magnetic 
moment, and magnetic anisotropy energy (bottom) of a single 
Fe adatom on Pt(lll) as a function of the number of self- 
consistent ly treated atomic shells around the adatom. 



have increased the number of the self-consistently calcu- 
lated perturbed atoms around the adatom from 12 (the 
first atomic shell around the impurity) up to 85 (the 
fourth atomic shell). It should be noted that the per- 
turbed atoms refer to substrate atoms and empty spheres 
(part of the vacuum region). 

In principle the magnetic anisotropy energy consists of 
two parts, namely the difference in total energy AE^ V 
and in the magnetic dipole-dipole interaction energy 
AEdd (shape anisotropy) between two given uniform ori- 
entations \i and v of the magnetization. Since for single 
adatoms the magnetic dipole-dipole interactions are of 
little importance only AE^ V is considered, which in turn 
was calculated by means of the force theorem as the cor- 
responding difference in band energies [15] , for details 
see Ref. [7 . As can be seen from Fig. [I] reliable con- 
vergence of the spin and orbital moments of the adatoms 
with respect to the number of shells of neighbors used 
was obtained if only two shells of atoms were taken into 
account. For the calculations of the magnetic anisotropy 
energy up to four shells of atoms were required to obtain 
converged values as illustrated in Fig. [I] 

The size of the magnetic moments of the adatoms ex- 
ceed those of the respective bulk materials Fe and Co, 
and of (complete) monolayers of Fe and Co deposited 
on Ir(lll) and Pt(lll). This behavior, now quite well- 
known to be characteristic for small magnetic clusters on 
top of metal substrates, [6], [4] is due to the lower coor- 




FIG. 2: Spin moments of the adatoms and induced spin mo- 
ments in the topmost substrate layers. The values correspond 
to the size of the moments along the easy magnetization axis, 
which is out-of-plane in all cases except for the Fe/Ir system. 

dination of the surface atoms which favors an incomplete 
quenching of orbital contributions. For both kinds of 
adatoms, Fe and Co, the magnetic moments are larger 
when deposited on a Pt substrate, a peculiar feature, 
which most likely is caused by the stronger polarization 
of Pt than that of an Ir substrate, see Fig. [2] In fact the 
induced spin magnetic moments in the nearest neighbour 
atoms of the Pt substrate are by about 0.02/i# higher 
than for the Ir substrate, and the polarization rapidly 
decreases by one order of magnitude for the second and 
third nearest neighbours. 

The spin and orbital moments of the adatoms are sum- 
marized in Table [I] At a first glance it can be realized 
that in relation to the bulk values the spin moments are 
considerably increased. Compared to the (theoretical) 
value of bulk Co (/^ cp = 1.6/xb) the spin moments of the 
Co adatom are increased by a factor of approximately 
1.3 if deposited on either substrate. For Fe, which has 
a bulk value of approximately /i^ cc = 2.1/iB, this ratio 
is with 1.6 slightly larger. While it is known that for 
3d bulk systems, LSDA density functional calculations 
predict the spin moments rather accurately (underesti- 
mating the experimental values by only about OT/ie), 
the orbital moments of Fe and Co, in particular, are off 
by about a factor of 2. Arguably, correlation effects may 
play a prominent role in predicting the size of the orbital 
polarization correctly [14 . The calculated values listed 
in Table [T] can therefore be expected to underestimate the 
actual size of the orbital moments of the adatoms. How- 
ever, it is worthwhile to consider the amount by which 
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FIG. 3: Contributions to the MAE of the adatoms and the 
substrate. From left to right and top to bottom: Co/Pt, 
Fe/Pt, Co/Ir, Fe/Ir 
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FIG. 4: Contributions of the substrate atoms in the surface 
layer (S) and the two layers below (S-l, S-2) to the total MAE. 

these values are increased compared to the bulk values of 
LSDA calculations. The orbital moment of Co in bulk is 
0.078/iB and that of Fe 0.043/iB- Due to the reduced co- 
ordination and the different chemical environment these 
values are increased for both an Fe and a Co adatom on 
the Ir substrate by a factor of approximately 6.3. In the 
case of a Pt substrate the Co and Fe orbital moments are 
9.3 and 14.6 times larger, respectively. This increase of 
the orbital polarization is accompanied by an enhance- 
ment of the magnetocrystalline anisotropy. 

Within our method it is possible to calculate the con- 
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FIG. 5: Variation of the magnetic anisotropy energy of a sin- 
gle Co adatom on Ir(l 11) (triangles) and on Pt (111) (circles) 
as a function of the orientation of the magnetization with 
respect to z-axis, as specified by the polar angle G. 
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FIG. 6: Variation of the magnetic anisotropy energy 
of a single Fe adatom on Ir (111) (triangles) and on 
Pt (111) (diamonds) as a function of the orientation of the 
magnetization with respect to z-axis, as specified by the polar 
angle G. 

tributions of individual atoms to the total MAE. The val- 
ues obtained are illustrated in Fig. [3] for the adatom and 
atoms in the surface layer of the respective substrates, 
and hence we can evaluate that portion of the anisotropy 
energy which is attributed to the substrate. Interestingly 
we find a significant dependence of these contributions on 
the type of deposited adatom. If Fe is deposited on Pt 
the substrate contributes only about 5.5%, whereas in 
the case of an Co atom the Pt atoms add 11.9% to the 
total MAE. In contrast, the atoms of an Ir substrate con- 
tribute 15.3% for the Co adatom, and even 25.5% if the 
adatom is Fe. The latter case is the only instance when 
the preferred magnetization direction due to the MAE is 
perpendicular to the surface normal. In that case, inter- 
estingly, the major contribution of the substrate does not 
come from the atoms in the surface layer closest to the 
adatom, but from the subsurface layer (S-l, c.f. Fig. |3|. 

Independent of the substrate for a Co adatom AE^ 
predicts strongly an out-of-plane orientation of the mag- 
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adatom 


-2.655 


3.359 


3.341 


0.267 


0.243 




Co 


ML 


1.395 


1.893 


1.900 


0.126 


0.142 




adatom 


2.982 


2.008 


2.040 


0.427 


0.494 



TABLE I: Anisotropy energies [meV], spin and orbital mag- 
netic moments [/jLb] of single monolayers of Fe and Co 
on Pt(lll) and Ir(lll) as compared to the corresponding 
ad atoms values. 



netization with the easy axis along z direction, see Fig.[5j 
while for an Fe adatom, see Fig. [6j the choice of the sub- 
strate seems to be significant: on top of Pt(lll) an out- 
of-plane orientation with an easy axis along z applies, 
while deposited on Ir(lll) an easy axis along x is pre- 
ferred. 

Finally for matters of comparison to the single adatoms 
AE^ V of a complete Fe and a Co monolayer ferromagnet- 
ically coupled [15 to the substrate was studied within the 
framework of the fully relativistic SKKR method [7J . The 
results of this study are displayed in Table [T[ One immedi- 
ately observes that the orbital magnetic moments of the 
adatoms are more sensitive to the chemical environment 
than the spin moments and in the meantime, the orbital 
moment anisotropy is larger on the Pt substrate than on 
the Ir. In comparing now Figs. [5] and [6] with the values 
in Table |IJ one easily can see that AE^ U is substantially 
larger for a Co adatom on Ir(lll) or Pt(lll) than for the 
corresponding monolayer. In the case of Fe adatoms one 
even has a reversed situation: according to our calcula- 
tions and in a good agreement with experimental findings 
[T6] a single monolayer of Fe on Ir(lll) or Pt(lll) ex- 
hibits an in-plane magnetic anisotropy while a single Fe 
adatom on Pt(lll) is strongly perpendicularly oriented. 

An additional contribution to the magnetic 
anisotropy comes from the magnetic dipole-dipole 
interaction energies (shape anisotropy), which is 
always negative |17j and consequently favors an 
in-plane orientation of the magnetic moments. 
For single magnetic monolayers on metal sub- 
strates the shape anisotropy is rather small (typ- 
ically about -0.1 meV) and was not taken into 
account in the present calculations. It should be 
noted that the shape anisotropy becomes very im- 
portant indeed whenever the number of magnetic 
monolayers is increased and then very often is the 
cause for so-called reorientation transitions, see 
for example the discussion in Ref. |20j. The num- 
bers given in Table [I] compare the band energy 
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FIG. 7: "Spin-projected" density of states of a single adatom 
of Fe (shaded area, dashed line) and of a single Co adatom 
(black line) on Pt(lll) (top) and Ir(lll) (bottom). 



contribution to the magnetic anisotropy energy 
of a monolayer with that of a single adatom. 

Fig. [7] shows the total spin resolved density of 
states (DOS) for the Fe and Co adatoms on the 
Pt and Ir substrate, respectively. Note that the 
magnetic field was taken along the surface nor- 
mal in all cases. As can be expected the DOS 
appears to be very similar for both Fe and Co, 
with the spin-up DOS always completely filled. 
The spin-down DOS of Fe and Co are of almost 
identical shape, but the peak is shifted slightly 
to lower energies for Co to accomodate the ad- 
ditional electron. Comparing the DOS on the Pt 
and Ir substrates, one can identify a broadening of 
the DOS on Ir which results in the slightly lower 
spin moments, S z , of Fe and Co (c.f. Table [T]). 

As has been shown in Refs. |18| and [19J the 
increase in the orbital moments, as compared 
to the bulk or the monoloayer cases, is caused 
by the difference in the filling of the d a (a = 
xy, xz, z 2 , yz, x 2 — y 2 ) orbitals. 



It has been noticed quite some time ago [2TJ [22] that 
in interfaces of Fe films with layers of the 5d elements 
the induced orbital moments can violate Hund's third 
rule. As a consequence of this rule the spin and orbital 
moments should align antiparallel for a less than half 
filled shell and parallel for a shell more than half filled. 
Even though strictly valid only for atoms it seems that 
Hund's rules are also applicable to solids, however, with 
exceptions. Since both Pt and Ir possess a more than 
half filled d-shell, J = L + S has to be expected. The 
relative orientation of spin and orbital moments are ex- 
plored in Fig. [8] where the atoms in the surface layer that 
occupy sites in the vicinity of the adatom are shown. A 
indicates parallel and a '-' antiparallel alignment of 
the moments. Note that the spin moments align parallel 
with the spin moments of the adatom, as also illustrated 
in Fig. [2] However, this is not strictly the case. Albeit 
not shown in this work, there are induced spin-moments 
in the subsurface layer which align antiparallel to the 
spin of the magnetic adatom. Fig. [8] shows, that in a Pt 
substrate the spin and orbital moment are parallel with 
Fe on top, and only a few moments are aligned oppo- 
site with the Co adatom on top. The situation is very 
much different for the Ir substrate, where many more 
atoms show antiparallel alignments. The arrangement of 
the atoms exhibiting this anomaly shows a symmetry ac- 
cording to the hexagonal 2D lattice in the case of the Co 
adatom. However, for the Fe adatom this arrangement 
does not have the same symmetric pattern, because the 
easy magnetization direction is in-plane along the x direc- 
tion. Small changes in the band filling have been 
shown to lead to such a behaviour [22j. 

Co nanostructures on Pt(lll) were already studied in 
the past in terms of thin films of Co n or (CoPt) n super- 
structures [23], n denoting the number of atomic layers 
or repetitions and in the form of finite chains of Co atoms 
[24| [25] . In Ref. [23] it was claimed that for a single layer 
of Co on Pt(lll) an in-plane orientation of the magneti- 
zation is preferred with a very small anisotropy energy. 
In the spin dynamics study of Ref. [24], which is based 
on the same computational approach as used in here not 
only the value of the magnetic anisotropy energy agreed 
very well with experiment [1 , but also the direction of 
the canted magnetization. A recent study [26] of the 
structure of a thin film Co/Pti3/Co revealed that an fcc- 
type stacking of Co was more favorable than an hep-type 
stacking, the interlayer distance between the Co and the 
first Pt layer being reduced by 10.1% as compared to that 
of bulk Pt. Although a film with two magnetic surfaces 
cannot be compared directly with a semi-infinite sub- 
strate coated with a monolayer of a magnetic metal (in 
the case of a semi-infinite system the Fermi energy is al- 
ways that of the substrate, i.e., differs from that of a thin 
film), these results indicate that layer relaxation might 
in special cases be important for investigating magnetic 
anisotropics. 
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FIG. 8: Alignment of the spin and orbital moments in the sur- 
face layer of the Pt and Ir substrate. The sign '+' means that 
the induced moments are aligned parallel and '-' that they are 
aligned antiparallel. The central, black, atom is the magnetic 
adatom. From left to right and top to bottom: Co/Pt, Fe/Pt, 
Co/Ir, Fe/Ir 



Experimentally a study of Co nanostructures on 
Pt(lll) [26j [27] seems to lead to a rather complicated 
situation. In Ref. [26 an out-of-plane magnetization of 
Co wires and islands is found and - in order to explain the 
measured domain wall width in the wires - an effective 
anisotropy constant between 0.08 and 0.17meV/atom for 
atoms within an island is proposed. In Ref. [27] the 
edge atoms of small islands are made responsible for their 
uniaxial out-of-plane magnetization. Of course none of 
these experimental results can be compared directly with 
results for a smooth Co monolayer on Pt(lll) exhibit- 
ing two-dimensional translational symmetry. Finite Co 
nanostructures (single adatoms or finite wires) on Ir(lll) 
and on Pt(lll), however, do show a strong perpendicular 
anisotropy and support the experimental findings. 

Up-to-now the case of using Ir(lll) as substrate is 
much less studied. From a spin-polarized STM study 
of (incomplete) monolayers of Fe on Ir(lll) it was found 
[28] that the stacking of the Fe atoms seemed to be of 
fec-type. Furthermore, by applying an external magnetic 
field to the tip (and of course also to the sample) from 
the impact of this field Bergmann et al. [28] concluded 
that "the observed superstructure is of magnetic origin 
with an out-of-plane magnetization". In order to check 
whether or not an externally applied field can be the 
cause for a possible misinterpretation of experimental re- 
sults, the effect of such a field on the band energy part of 
the magnetic anisotropy energy was simulated by consid- 
ering (non-selfconsistently) a Kohn-Sham-Dirac Hamilto- 
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FIG. 9: Variation of the band energy part of the magnetic 
anisotropy energy of a single monolayer of Fe on Ir(lll) as 
function of an applied (constant) external magnetic field. 

nian with the total magnetization pointing (a) along the 
surface normal, 

H(r) = col ■ p+f3mc 2 +V eff (r)I A +f3Z z [B e Jf(r) + B ext ] , 

and (b) in-plane. In the above equation a and (3 are 
Dirac matrices, J2 Z is the ^-component of the so-called 
spin operator, I4 is a four-dimensional unit matrix, and 
B ext is a small constant external (magnetic) field. This 
simulation is displayed in Fig. [9] and very clearly shows 
that the fields applied in experiment most likely only 
of marginally affect the size of the magnetic anisotropy. 
Possible sources of discrepancies between the weak 
in-plane anisotropy found in our calculations and the 
out-of-plane magnetization reported in [28] can be of 
quite different origin. For a smooth, two-dimensional 
translationally invariant Fe overlayer on Ir(lll) the 
anisotropy energy is only slightly negative. Therefore 
layer relaxations as well as the finite size of the sample 
might be of importance. Complicated geometrical dis- 
tortions due to the incompleteness of the atomic layers 
are very difficult to take into account theoretically. Yet 
another possibility is that the experimentally found 
out-of-plane magnetization in Fe/Ir(lll) is not caused 
by a perpendicular anisotropy but rather by other 
factors like, e.g., a complicated chiral rotation of the 
magnetization due to higher order exchange interactions 
[29] . It might even turn out that perfect monolayers 
and single adatoms such as considered in here cannot 
reflect sufficiently well the actual situation mapped in a 
particular experiment. 
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